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Abstract 

We review some pseudo-planar deformations for the equations of 
incompressible isotropic nonlinear elasticity first introduced in 1985 
by Rajagopal and Wineman. We extend this class of deformations 
to compressible isotropic and transverse isotropic materials, and also 
consider the influence of gravity. We consider some new approximate 
solutions and we discuss the possible relevance of such solutions to 
the understanding of the complex structure of the fields equations of 
nonlinear elasticity, using weakly nonlinear theories. 
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1 Introduction 



The starting point of our discussion is the 1985 paper by Rajagopal and 
Wineman entitled New Exact Solutions in Non-Linear Elasticity. In that 
paper some new exact solutions to boundary value problems of nonlinear 
elasticity were established. These solutions constitute an important break- 
through point in the literature and to appreciate this fact, we have to review 
what was going around in Continuum Mechanics and in the field of nonlinear 
elasticity at the time. 

The origins of the modern theory of nonlinear elasticity must be related 
to the pioneering works of Rivlin, Green, and Adkins just after World War II. 
These fellows were the first to elaborate a reasonable and general notation to 
derive the right balance equations for nonlinear elasticity, and were the first 
to determine classes of exact solutions in nonlinear elasticity (for a summary 
of this earlier work see Green and Adkins, 1970 fl 

Then in 1955, Ericksen proved that homogeneous deformations are the 
only controllable static deformations possible in every hyperelastic material. 
A controllable deformation is a deformation that is produced in a material 
by the application of surface tractions alone. A controllable deformation 
that can be effected in every homogeneous isotropic material is referred to 
as a universal solution. Ericksen's result concerning universal deformations 
has had a profound influence on the development of nonlinear elasticity. For 
many years afterwards, there was "the false impression that the only defor- 
mations possible in an elastic body are the universal deformations" (Currie 
and Hayes, 1981). This is not exactly true and for a list of notable excep- 
tions see for example the book by Ogden (1984). Currie and Hayes (1981) 
were right to point out that after some initial interest in the search for exact 
solutions in nonlinear elasticity, there followed a long period of inactivity 
concerning this enterprise. The paper by Currie and Hayes (1981) was the 
seed necessary to revitalize the search for possible solutions beyond universal 
solutions. The 1984 and 1985 papers by Rajagopal and Wineman were the 
first relevant new fruits. 

In 1979, Ratip Berker proposed some new exact solutions for the Navier- 
Stokes equations. These are an infinite set of nontrivial solutions for an in- 
compressible viscous fluid contained between the two parallel infinite plates 

1 To be more precise, Gabriella Armanni published a note in 1915, transmitted by Vito 
Volterra to the Nuovo Cimento journal [9 (1915) 427-447], with some radial solutions for 
a spherical nonlinear elastic solids, but this was an isolated episode. 
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rotating with constant angular velocity around the fixed normal axis. These 
flows are pseudo-plane flows of the first kind (Berker, 1963). Rajagopal was 
able to extend these solutions to a huge class of non-Newtonian fluids (Ra- 
jagopal, 1992) in an elegant note (Rajagopal, 1982). Rajagopal and Wineman 
(1984, 1985) considered the solid mechanics counterpart of the Berker solu- 
tions to obtain exact non-universal solutions for all incompressible isotropic 
elastic materials. This family of deformations is then generalized in (Ra- 
jagopal and Wineman, 1985) by considering that the bottom and top plates 
rotate with different angular velocities. These deformations are admissible 
solutions for all Mooney-Rivlin incompressible elastic materials. 

Rajagopal has been a fine reader and connoisseur of the Berker encyclo- 
pedic article (Berker, 1963) where the pecularities of pseudo-plane flows are 
discussed in great detail. Moreover, he established the formal analogy be- 
tween the Navier-Stokes equations and the equations determining admissible 
deformations for neo-Hookean solids, and thus turned the Berker article into 
a formidable source of inspiration for new classes of solutions in nonlinear 
elasticity. 

To illustrate why pseudo-plane deformations constitute an interesting 
class of potential solutions in nonlinear elasticity, we consider some quan- 
titative details. With material and spatial rectangular Cartesian coordinates 
(X,Y,Z) and (x,y,z), respectively, a pseudo-plane deformation is a defor- 
mation of the form 

x = x(X,Y,Z), y = y(X,Y,Z), z = Z. (1) 

An interesting way to generate such a deformation is to consider a plane 
deformation solution which depends on a certain number of arbitrary con- 
stants, a,b,c, . . ., say, and to replace these constants with functions of Z. For 
example, by going from 

x = x(X,Y;a,b,c), y = y(X,Y;a,b,c), z = Z, (2) 

to 

x = x(X, Y; a{Z), b(Z), c{Z)), y = y(X, Y\ a{Z), b(Z), c(Z)), z = Z. 

(3) 

This idea was independently considered also by Hill and Shield in 1986 but 
only for neo-Hookean materials. It is in Rajagopal's work that it has resulted 
in many new exact solutions of non-linear elastic layers. The Berker solution 
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is reconsidered also in (Fosdick et al, 1986); a non-uniform extension of a 
slab of Mooney-Rivlin material is considered in (Rajagopal et al., 1986); 
the rectilinear deformation of a general incompressible slab is considered in 
(McLeod et al., 1988); and not only deformations but also motions of elastic 
slabs are considered in (Carroll and Rajagopal, 1986). For pseudo-plane 
deformations and motions of the second kind (Berker, 1963), we refer to 
(Hayes and Rajagopal, 1992) and (Horgan and Saccomandi, 2003). In the 
paper (Rajagopal, 1984), the Berker (1979) solution is generalized to the 
case where the Newtonian flow is contained between the two parallel infinite 
porous plates. This suggests a generalization of ([TJ such that 

x = x{X,Y,Z), y = y(X,Y,Z), z = \(Z). (4) 

This class of deformation is used in (Saccomandi, 2005) to obtain new exact 
solution for neo-Hookean solids. 

Pseudo-planar solutions of the Navier-Stokes equations have not been 
the only inspiring source of analogies to obtain new solutions in nonlinear 
elasticity from existing solutions in fluid mechanics. By considering the cel- 
ebrated Jeffrey and Hamel convergent and divergent flows in intersecting 
planes (Berker, 1963), Rajagopal and co-workers have produced interesting 
deformations admissible in nonlinear elasticity (Fu et al., 1990; Tao and 
Rajagopal, 1990; Rajagopal and Carroll, 1992; Rajagopal and Tao, 1992; 
McLeod and Rajagopal, 1999). This class of deformations has also been 
considered by Klingbiel and Shield (1962), but Rajagopal has been able to 
reverse the coin of Continuum Mechanics to transfer his fluid mechanics ex- 
pertise to solid mechanics. 

The aim of this note is to push a little further the ideas of Rajagopal and 
investigate what happens to pseudo-planar solutions in the compressible case 
of nonlinear isotropic and anisotropic elasticity. Hence in Section 13. 1[ we use 
(j3J) to study rectilinear shear deformations coupled to an axial deformation 
of an elastic slab in the framework of compressible materials. We show that 
if we take into account the weight of the elastic slab, then the rectilinear 
shear solutions are non-trivial. It is usual in nonlinear elasticity theory to 
assume that the weight of an elastic structure is negligible when compared 
to elastic forces, but here we find that for certain fields of application such 
as geophysics, this is not necessarily the case, especially for thick layers. 
Then for incompressible materials (Section l3.2p . we derive a universal relation 
that went unnoticed by McLeod et al. (1988) and which proves very useful 
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in simplifying the analysis of the problem (For universal relations in the 
framework of rectilinear shear see (Saccomandi, 1996)). We use this class of 
deformations further in Section H] to generalize some results to the case of the 
rectilinear shear of a transversely isotropic elastic material. In (Destrade et 
al., 2009) we considered orthotropic materials were the fibers are arranged 
in a special plane; here, the fibers may have a general orientation in space. 

In order to make progress and obtain exact solutions in nonlinear elas- 
ticity, a special constitutive assumption must be enforced. Indeed only Er- 
icksen's universal solutions are valid irrespective of the choice of the strain 
energy density function. Sometimes the constitutive assumptions that have 
been considered in the literature were followed only for reasons of mathemati- 
cal convenience. This approach has generated a huge number of strain energy 
density functions that have no real connection with experimental evidence. 
Here we use a different approach, using the so called wealkly nonlinear theory 
of elasticity. The strain energy associated with the classical linear theory of 
elasticity is of second order in the Green strain; here to investigate the non- 
linear behavior of elastic materials, we use strain energies that are of the 
third and fourth order in the strain. In such a way we are able to bring out 
some interesting general features of nonlinear elastic materials. For example, 
we show in Section H] that the coupling among the various mode of deforma- 
tions occurs at lower order of non-linearity for anisotropic materials than for 
isotropic materials. 



2 Governing equations 

We call x(X) the current position of a particle which was located at X in 
the reference configuration. Two kinematic quantities associated with this 
motion are the deformation gradient and the left Cauchy-Green strain tensor, 

F = dx/dX, B = FF T , (5) 

respectively. In Section [3X we consider hyperelastic, isotropic materials, and 
so we introduce the strain energy density W = W(Ii, I2, 13), where I\ and I2 
are the first and second principal invariants of B, respectively, given by 

I, = tr B, l 2 = \ [If - tr(B 2 )} , (6) 
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and I 3 = J 2 where J = det F. The general representation formula for the 
Cauchy stress tensor T reads 

T = (3 I + f3 1 B + (3„ 1 B- 1 , (7) 

where 

/3 = 2 J" 1 [I2W2 + I3W3] , /3 1 =2J" 1 iy 1 , /3„ 1 = -2Jiy 2 , (8) 

and Wi = dW/dh, W 2 = dW/dlz, W 3 = dW/dI 3 . 

If the material is incompressible then the only admissible deformations 
are isochoric, i.e. J = 1 at all times, so that I 3 = 1, and W = W(I\, I 2 ) only. 
In this case, 

T = -pi + 2W 1 B - 2W 2 B~\ (9) 

where p is the yet indeterminate Lagrange multiplier introduced by the con- 
straint of incompressibility. 

The link with linear elasticity is made by defining /x, the infinitesimal 
shear modulus, as 

fi = 2(W 1 + W 2 )\ Il=l2=3j3=1 . (10) 

In Section HJ we consider a special class of hyperelastic anisotropic mate- 
rials, namely transversely anisotropic materials, i.e. with a single preferred 
direction. We call M a unit vector along this direction in the reference 
configuration and we introduce the anisotropic invariant 

I 4 = M ■ CM- 1 = m ■ m - 1, (11) 

where C = F T F is the right Cauchy-Green strain tensor, and m = FM. 

We restrict our attention to compressible, transversely isotropic materi- 
als for which W = W(Ii, I 2 , h, h)- Then the constitutive equation ([7]) is 
replaced by 

T = (3 I + PiB + P^iB" 1 + /3 4 m <g> m, (12) 

where /?4 = 2J~ 1 W^ W4 = dW/dI±. For incompressible, transversely isotropic 
materials with W = W(h, I 2 , J4), the representation formula ([9]) is replaced 
with 

T = —pi + 2W\B — 2W 2 B~ 1 + 2W4T71 ® m. (13) 

The balance equation of linear momentum, in the absence of body forces, 
reads 

divT = 0. (14) 
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If we consider the presence of body forces such as gravitjo, for example, then 
we have 



where pg is the weight per unit of volume of the body in the current config- 
uration and p is the current mass density. By the conservation of mass, 



where po is the mass density in the reference configuration. 

3 Rectilinear shear and axial stretch 

Among the class of deformations in a special status is detained by the 
deformations 



x = AX + BY + f(Z), y — CX + DY + g(Z), z = h(Z), (17) 



where A, B, C, D are constants and /, g, h are functions of Z alone. These 
deformations consist of two shearing deformations f(Z) and g(Z) in the Z- 
direction, combined to an inhomogeneous stretch h(Z) along the Z-axis, and 
superimposed on a plane homogeneous deformation. 

The class of deformations (IP7|) is not universal, but it reduces the balance 
equations to a system of three ordinary differential equations for any choice 
of the strain energy density. Moreover, at A = D — 1, B — C — 0, it is 
the static counterpart to the usual longitudinal and transverse wave motions 
propagating along the Z-axis and polarized along the three reference axes. 
For the choices h(Z) = \Z, A = D = A" 1/2 , B = C = 0, where A is a 
constant, we obtain an isochoric deformation first considered in a dynamical 
context by Carroll (1967), and then by Rajagopal and co-workers. In those 
papers, the material was incompressible, while here we are considering both 
compressible and incompressible materials. Another interesting case occurs 
for the choices 



A = D = cosn, £ = -sinfi, C = sinfi, h(Z) = XZ, (18) 



2 In nonlinear elasticity we do not know of exact solutions taking into account the 
presence of gravity. In linear elasticity and in incremental elasticity, a classical exact 
solution taking gravity into account is the stretching of a bar by its own weight, see the 
book by Timoshenko and Goodier (1982) and the article by Fosdick and Shield (1963), 
respectively. 



divT + p# = 0, 



(15) 



Po = Jp, 



(16) 
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where f2 and A are constants. Here we recover the field studied by Rajagopal 
and Wineman (1984, 1985), describing the deformation of a nonlinear elas- 
tic layer contained between two infinite parallel rigid plates, each of which 
undergoes the same finite rotation. 

Considering now the class of deformations f|T7|) in all generality, we find 



A B f 
C D g' 
h! 



[F- v 



J 



-i 



Dh' -Bti Bg' - Df' 
-Cti Ah' Cf - Ag' 
AD-BC 



(19) 

where J = (AD — BC)h' and the prime denotes differentiation with respect 
to Z. Also, 



\B] 



A 2 + B 2 + f' 2 
BD 
fh' 



AC + BD + f'g' C 2 + D 2 + g' 2 



AC + BD + f'g' fh' 
g'h' 

g'h' h' 2 



(20) 



and 



31, 



[B- 1 ]^ J-^ + D 2 )/!' 2 , 

[B-% = -J~ 2 (AC + BD)h' 2 = [B-%, 
[B-% = r 2 [(AC + BD)g' - (C 2 + D 2 )f']h' = [B 
[B- l ] 22 = r 2 (A 2 + B 2 )h' 2 , 

[B-%3 = ,r 2 [(AC + BD)f - (A 2 + B 2 )g']h' = [B-% 2 , 
[B-% 3 = J- 2 [{Bg' - Df) 2 + {Ag' - Cf') 2 + (AD - BCf 



(21) 



The principal invariants follow as 

h = A 2 + B 2 + C 2 + D 2 + f 2 + g' 2 + h' 2 , 
h = (AD - BC) 2 + (Df - Bg') 2 + (Ag' - Cf') 2 
J 3 = (AD - BC) 2 h' 2 . 

Now we may compute the Cauchy stress components for compressible 



(A 2 + B 2 + C 2 



D 2 )h' 2 , 
(22) 
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solids from ([7]) as 



T u = p + /3i(A 2 + B 2 + f' 2 ) + .r 2 {C 2 + D 2 )h' 2 , 

T 22 = (3 + (3i(C 2 + D 2 + g' 2 ) + /?_! J~ 2 {A 2 + £? 2 )/i' 2 , 

T 33 = /3o + M' 2 + J- 2 [(^' - -D/') 2 + (V - Cf) 2 + (AD - 5C) 2 ], 

T 12 = ft (AC + + /V) - f3^J~ 2 l(AC + BD)h' 2 , 

Tia = Pifti + J- 2 [(AC + BD)g' - (C 2 + D 2 )f]ti, 

T 23 = Pxg'ti + J~ 2 [(AC + B£>)/' - (A 2 + B 2 )g']h'. (23) 

For incompressible solids, we may use the same formulas, by replacing /3 
with —p, (3\ with Wi, and j3-i with W 2 - 

3.1 The compressible case 

Now we write down the balance equations in the absence of body forces, 
equations (fT4"j) . They reduce to 

T[ 3 = 0, 7*3 = 0, T 3 ' 3 = 0, (24) 

and they may readily by integrated as 

Ti 3 = ki, T 23 = k 2 , T 33 = k 3 , (25) 

where fci, fc 2 , k 3 are constants. This is a non-homogeneous quadratic system 
of three equations with constant coefficients for the three unknowns g', h! . 
Once solved, it yields /', g', h' as constants, showing that the only possible 
solutions are the trivial homogeneous solutions. This is in contrast with the 
dynamic counterpart of these deformations (see Destrade and Saccomandi, 
2006). 

Next, consider the case where the gravity body force is present, so that 
pg has component — pg = —p$J~ l g along the Z-direction. Then the balance 
equations (TT5T) reduce to 

7^ = 0, 7^ = 0, Ti^poiAD-BCy'g, (26) 

and they may readily by integrated as 

T 13 = k u T 23 = ~k 2 , T 33 = p (AD-BC)- 1 gZ + k 3 , (27) 
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where k\, k 2 , k^ are constants. In contrast to (I25p . this differential system 
admits non-homogeneous solutions. 

Take the case A — D — 1, B = C = 0. Then, the governing equations 
fl27|) reduce greatly, to 

[^-^x{h')- 2 ]fh' = k u 
\Pi-^i{h'r 2 }g'h' =~k 2 , 

0o + m') 2 + /M/0~ 2 [i + iff + (</) 2 ] = PogZ + ~h- (28) 

To set down a boundary value problem, we consider a slab of finite thick- 
ness H along Z and of infinite extent otherwise. Then we scale the lengths 
with respect to H and the /3's with respect to fi, the infinitesimal shear mod- 
ulus defined in (TTUj) . The non-dimensional version of the system (I28p is thus 
written in the domain ^ Z ^ 1, as 

Qif' = h, Qig' = k 2 , (Q a + Q 2 )ti = (p gH/n)Z + k 3 , (29) 

where 

Q, = 2(Wi + W 2 )/fi, Q 2 = 2(W 2 + W 3 )/n, h = k/fi. (30) 

Notice that when /' ^ 0, g' ^ 0, we can deduce from fl29|) i p that /' is 
proportional to g', thus reducing the dimension of the system of ordinary 
differential equations. Notice also that we can differentiate fT251 with respect 
to Z, to get the equivalent form of the governing equations, 

(Qif)' = 0, (Qig')' = 0, [(Qi + Q2)h']' = ( Po gH/fi). (31) 

It is possible to impose Dirichlet boundary conditions by prescribing the 
displacements at the bottom and the top of the slab, or Neumann boundary 
conditions by imposing the values of shear stresses T\ 2 , T13 and of the normal 
stress T 33 on the faces of the slab. Mixed boundary conditions are also 
possible. 

Now we specialize the analysis to the Murnaghan strain-energy density: 

A + 2u t9 u t / + 2m o m T T n T . , 

W = -J^Jl + f J2 + -34- ^1 + 7^1^ + g J 3 , (32) 

where A and fi are the Lame coefficients of second-order elasticity, and /, m, n 
are the Murnaghan third-order constants. Here, Ji, J 2 , J3 are another set of 
independent invariants, related to the principal invariants of strain through: 

J 1 = h-3, J 2 = 2/ 1 -/ 2 -3, J 3 = J3-J2 + /1-I. (33) 
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For this strain energy, 



„ A + 2/i + m T / + 2m r2 m T 

„ _ A + 2/i l + 2m 2 m 

<3l + * = — Jl + — J? + ^ 2 - (34) 

To work out an explicit example, we take the deformation to be in the 
form f(Z) = u(Z), g(Z) = 0, h(Z) = Z + w(Z), for some functions u, w of 
Z, so that the mechanical displacement is x — X = [u, 0, wf. Then we find 
that 

J 1 = 2w' + (u') 2 + (w') 2 , J 2 = (u') 2 , J 3 = 0. (35) 

We also impose the displacements boundary conditions /(0) = 0, h(0) = 1, 
h(l) = £, equivalent to 

n(0)=0, w(0)=0, w(l) = £-l. (36) 

First of all we make the connection with the solution of linear elasticity. 
There, u = u , w = w , say, where u and w are infinitesimal quantities, 
in the sense that Mq, and higher powers are negligible when compared to 
|wo| an d \wq\. At that order of approximation, (131!) reduces to 

u = 0, w = — , (37) 

with solution: 

Uo = kiZ t Wo =P£l!—(Z-l)Z+(£-l)Z, (38) 

where the integration constant k\ is to be determined later on, from the 
the stress boundary condition imposed on T 13 (l). It is clear that there is a 
limitation to this solution in the sense that the expression for wq must stay 
within the so-called "elastic limit" , and this is not guaranteed for slabs with 
large thickness H. In real world applications where gravity plays a significant 
role (such as geophysics), it could well be the case that a non-linear correction 
is required to account for the effects of a large layer thickness. 

We now move on to the next order of elasticity theory, that is the one 
encompassed by the strain energy f )32|) . To obtain a solution for the bound- 
ary value problem, we use a simple perturbation method, by expanding the 



11 



displacement as u = u + U\ + . . ., w = w + w x + . . ., say, in the same spirit 
as the method of successive approximations introduced by Signorini (1949), 
and later studied by many authors (see for example, Lindsay (1985)). Here, 
«0) are given by fl38|) . and \w\\ are of the same order as and Wq, 
and higher orders are neglected. From (l3Tj) we obtain the following system 
of equations, 



A + 2/i + m . . 

-u w Q 



0. 



w[ + 



A + 2/i + m / 2 3A + 6/i + 2/ + 4m 



2p 



w 



/\2 



0. 



(39) 



Clearly now, each component of the next order solution involves a combina- 
tion of the longitudinal deformation (through wq) and the shear deformation 
(through uq). Owing to the form ( |38l) of the lower order solutions, these 
equations reduce to 



u' + 



w>{ 



A + 2/i 



m 



/i(A + 2//) 



hpgH = 0, 



3A + 6/i + 21 + 4m 
2/i(A + 2/i) 99 



pgH 



_2(A + 2/x) 



(2Z 



0. (40) 



Hence, the rectilinear shear is no longer a simple shear deformation. Here 
we have, according to ( T40l) i . a quadratic variation of u with Z, 



u 



kiZ 



A + 2/i + m pgH 
A + 2/i 2/i 



(41) 



and, according to fT4TJ]) 9. a cubic variation of w with Z, which we do not 
reproduce for brevity. 

In his treatment of gravity seismic waves Biot (1940) gives pg/p as being 
of the order of 2 x 10~ 6 m -1 for the Earth. This means that the non-linear 
correction above is of the same order as the linear solution when H ~ 1000 
km, or about one-sixth of the Earth radius. The value chosen by Biot is 
typical of common rocks such as granite, but it is unlikely to be representative 
of material properties at such depths. For softer grounds such as mud or 
sediments, pg/p is much larger, of the order of 2.5 x 10~ 4 m _1 say (Holzer 
et al., 2005), indicating that H ~ 8 km for the non-linear correction to be of 
the same order as the linear solution. At such depth however, consolidation 
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has occurred and the medium is much stiffer than near the surface. The 
conclusion of these estimates is that on physical grounds, we do not need 
to push the Murnaghan expansion to fourth order, because third-order non- 
linear effects are simply a very small correction to the linear solution. 

Nonetheless, simple shear deformations are an important component of 
most geophysical applications. For e xample, they form the bas is of an ex- 
planation of the folding phenomenon (IManz and Wickhaml . Il979l ). Similarly, 
the determination of shear strength, viscosity, and internal friction data for 
deep crust and upper mantle rock and mineral analogues under geophysically 
realistic conditions of very high temperatures and pressures is required in or- 
der to interpret eart h quake origins, seismic signal generation, and explosions 
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Okumura et al.l (120101 ) and references therein) a s well as the rheolog ical 



evolution of the microstructure of mantle materials (IKarato et al.l . Il998l ). It 



is clear that in order to have a realistic geophysical picture, it is necessary 
to conduct a numerical investigation of equations fl3Tl) . and to include the 
effects of very high temperatures and pressures. We note that under these 
circumstances, the constitutive prope rties of geo physical materials should 
quite softer than those indicated by iBiot fll940h . leading to a critical H 
much more smaller than the one we obtained above in our crude estimation. 
This means that the effect of gravity, and the corresponding inhomogeneous 
correction predicted by the Murnaghan theory, might both have to be taken 
into account after all. 



3.2 The incompressible case 

Here we go back to the incompressible already examined by 

McLeod, Rajagopal, and Wineman (1988). The aim of this subsection is to 
show that the use of a universal relation reduces the problem to a simple 
formulation. 

In the incompressible case the deformation (TTT|) must be isochoric. Here, 
we follow the choice of Carroll (1967) and of Mc Leod et al. (1988), by 
taking h(Z) = XZ, A = D = A~ 1/2 , B = C = 0. We study the equations of 
equilibrium in the absence of body forces. In contrast with the compressible 
case, non-trivial (non-homogeneous) solutions are possible, due to the effects 
of the Lagrange multiplier p. Indeed, the first two equations of equilibrium 
lead to dp/dx = ci, dp/dy = C2, where c\ and c 2 are constants. Assuming 
that Cic 2 7^ excludes the possibility of trivial homogeneous deformations. 
Then, the determining equations for the rectilinear shear deformation f(Z) 
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and g(Z) are 

Q x f = Cl Z + k u Q ig ' = c 2 Z + k 2 , (42) 

while the third equation, corresponding to fl26l) . determines the unknown La- 
grange multiplier p. We recall that gravity is a conservative external force: 
it would change the form of p but would not influence the shearing deforma- 
tions, which is why we did not include it in this section. 

We note that McLeod et al. (1988) did not make use of the following 
universal relation 

obtained directly from (I42p . With this universal relation, the discussion of 
boundary value problems is greatly simplified, as we now see for the most 
general case of fourth-order incompressible elasticity, 

W = \i tr (E 2 ) + - tr (E 3 ) + V (tr(£ 2 )) 2 , (44) 
3 

where E = (F T F — I)/2 is the Green strain tensor, and /i, A, and V are 
second, third-, and fourth-order elasticity constants, respectively (Hamilton 
et al., 2004; Ogden, 1974). For convenience, we work with the equivalent 
strain energy density 

W(h, h) = C 10 (/i - 3) + C 01 (I 2 - 3) + Cao(/i - 3) 2 , (45) 

where Ci , C i, and C 20 are constants. At the same degree of approximation 
in the Green strain E, it covers (jHJ) with the identifications (Destrade et al., 
2010): 

fji = 2{C w + C 01 ), A = -8(C W + 2C 01 ), V = 2(C 10 + 3C 01 + 2C 20 ). 

(46) 

and at C 2 q = 0, it also covers the Mooney-Rivlin case of McLeod et al. 
(1988). Then the determining equations (l4"2"j) reduce to a coupled system for 
the shear modes: 

[a + P(f 2 + g' 2 )]f = c 1 Z + k 1 , 

[a + (3(f' 2 + g' 2 )}g' = c 2 Z + k 2 , (47) 
where the scalar a and (3 are defined as 

a = l + (A 2 + 2A- 1 -3)/3, = 2 °™ n . (48) 

^10 + ^01 
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We consider the following boundary conditions /(0) = /(l) = and 
g(0) = g(l) = 0. Using symmetry arguments (see Saccomandi (2004) for 
details), it may then be shown that k\ and k 2 are such that 

[a + /3(/' 2 + g' 2 )]f = ci (Z-i), 

[a + /?(f 2 + </ 2 )K = c 2 {Z-\). (49) 

Using now the universal relation, which reads /' = (ci/c 2 )g', the problem is 
reduced to that of a classical rectilinear shear deformation 

a(g') +p(l + ^\ {g'f = c 2 (Z-\), g(0) = g(l) = 0. (50) 

Note that at C 2 o = we recover the solution of McLeod et al. (1988) for the 
Mooney-Rivlin material: 

g(Z) = ^Z(Z-l). (51) 

Now we perform a perturbation scheme, as in the previous section. Hence 
we take g in the form g — Vq + V\ + v 2 + . . ., say, where vq is infinitesimal, 
v\ is of order Vq, v 2 is of order Vq, etc., and we take the stretch in the form 
A = l + e + e 2 + ..., where the elongation e is infinitesimal. Then we find 
that the quantity a expands as a = 1 + 3/3e 2 + . . . and /3 remains the same. 
We also find that the single ordinary differential equation ( |50|) gives in turn, 

v' = c 2 (Z-l), so that v {Z) = ^Z{Z - 1), (52) 

then 

v[ = 0, so that V\ = 0, (53) 

showing that we were right to push the expansion of W to the fourth-order, 
and finally, the non-linear correction, 



(l + |) K) 3 + 3e 2 K) 



(54) 



We thus find the solution for g as 



g(Z) = °jZ(Z - 1) jl - I [(cl + cl)(2Z 2 -2Z+1) + 12e 2 ] +...}, (55) 



and for /, by integrating the universal relation, as f(Z) = (cx/c 2 )g(Z). 
Thanks to the successive approximation solutions, and to the universal rela- 
tion, we have been able to understand the role of higher order non-linearity 
and the coupling between the two shear components of the pseudo-planar 
deformation. 
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4 Rectilinear shear in anisotropic materials. 

Let us consider a material reinforced with fibers aligned with the direction 
M = M 1 E 1 + M 2 E 2 + M 3 E 3 , where the constants M u M 2 , M 3 are of the 
form 

Mi = sin 6 sin <ft, M 2 = sin# cos0, M 3 = cos8, (56) 

with 9, the elevation and 0, the azimuthal angle. 

We restrict our attention to the incompressible case and the isochoric 
motions 

x = X + f(Z), y = Y + g(Z), z = Z. (57) 
Then m = FM has components 

[m]i = [Mi + f'M 3 , M 2 + g'M 3 , M 3 f , (58) 

so that 

h = m • m - 1 = [(f) 2 + {g'f] M 3 2 + 2(f'M 1 + g'M 2 )M 3 . (59) 

The representation formula f lTS]) gives the following shear stress compo- 
nents 

T 13 = 2{W X + W 2 )f + 2W 4 M 3 (M 1 + f'M 3 ), 

T 23 = 2(Wi + W 2 )g' + 2W 4 M 3 (M 2 + g'M 3 ). (60) 

In passing we note that although W may be linear in I\ and I 2 (as in the 
Mooney-Rivlin model), it may not be linear in because otherwise there 
would be non-zero shear stresses in the reference configuration. The balance 
equations (TH|) reduce now to 

-p x + T[ 3 (Z)=0, - Py + T^(Z) = 0, n 3 (Z) = 0. (61) 

The third equation here can be solved by an appropriate choice for p. The 
first and second equations are compatible when 

T[, = 2c x , ^3 = 2c 2 , (62) 

where c% and c 2 are constants. Integrating, we obtain 

(W 1 + W 2 )f + iy 4 M 3 (M 1 + f M 3 ) = ciZ + c 3 , 

(W 1 + W 2 )g' + W 4 M 3 (M 2 + g'M 3 ) = c 2 Z + c 4 , (63) 
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where c 3 , c 4 are constants. 

First, consider the case of a single shearing deformation: /' ^ 0, g' = 
say. Then the second equation above results in 



W 4 M 3 M 2 = c 2 Z + c 4 , (64) 

which cannot be satisfied in general because W4 is not constant. Therefore 
a shearing deformation in a single direction is not compatible with a general 
orientation of the fiber arrangement. Next, consider the combination of two 
shearing deformations f'g' 7^ 0. Clearly then, the governing equations (1631) 
are compatible with a general fiber distribution in space. This observation 
also holds when the angles 9 and in (|56|) depend on Z: 9 = 9(Z) and 
= 4>{Z). This has of course important biomechanical implications, because 
it is well established that collagen fiber bundles change their orientation in 
soft tissues (Holzapfel and Ogden, 2009). 

Now consider the standard reinforcing model for the Mooney-Rivlin solid, 

W = C 10 (/i - 3) + C m {I 2 - 3) + i T (/ 4 - l) 2 , (65) 

where C\q, Cqi, 7 are constants. Then the governing equations (1631) read 



Vf' + 7[(/' 2 + 9' 2 )Mt + (Mi/' + M^M^M,^ + /'M 3 ) = c x Z + c 3 , 
^' + 7[(f 2 + ^' 2 )M 3 2 + (M 1 f + M 2 g')M 3 ]M 3 (M 2 + g'M 3 ) = c 2 Z + c 4 , 

(66) 

where \i = 2(C W + C i). The isotropic version (at M 3 = 0) of this system is a 
non-homogeneous linear decoupled system of equations for /' and g', which 
may be solved for one shear deformation independently of the other. Clearly, 
the anisotropy gives a coupled, non-linear (cubic) system of equations for the 
shear deformations. We may say that anisotropic materials are more sensitive 
than isotropic materials to nonlinear effects. For a detailed inve s tigati on of 



this type of sensitivit y, we refer to the works of iMerodio et al.l (120071 ) and 



Destrade et all (120091 ). 



5 Concluding remarks 

We revisited a class of deformations previously studied into details by Ra- 
jagopal and coworkers over several papers. This class of deformation is the 
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pseudo-planar version of some homogeneous deformations, composed of two 
inhomogeneous rectilinear shears superimposed onto an homogenous defor- 
mation. 

First, we extended the class of deformations from incompressible to com- 
pressible materials, by introducing the possibility an inhomogeneous axial 
stretch. In that case, the class of deformations is clearly not isochoric. In the 
presence of a body force such as gravity the solutions to the problem under 
investigation are not trivial i.e. we have non-homogeneous deformations. In 
the linearized case, the axial stretch and the rectilinear shear functions are 
uncoupled, but in the weakly nonlinear theories of elasticity, the situation 
becomes more intricate. At the order above linear elasticity (i.e. third-order 
elasticity), we saw that the rectilinear shears are coupled to the axial stretch, 
although the two rectilinear shears do not interact with each other. Presum- 
ably, we would find a fully coupled situation at the next order (fourth-order 
elasticity), but it was not necessary to go that far, because we computed that 
those effects are negligible for commonly used values of physical constants. 
This situation nonetheless revealed an important feature of the structure of 
the equations of nonlinear compressible elasticity. Indeed, volume variations 
are very important in the general nonlinear theory. Therefore a special math- 
ematical compressible model which supports isochoric deformations must be 
seen as the exception rather than the rule, and must be handled with care 
because it might give false informations about the physics of compressible 
elastic materials. 

In the incompressible case we noticed an important universal relation, 
which allowed us to simplify the qualitative analysis of the determining equa- 
tions for the rectilinear shear unknowns. We provide a simple approximate 
solution of this problem. 

Then we considered what happens in the framework of anisotropic mate- 
rials, where lower-order coupling appears between the shearing deformations. 
This coupling is due to the presence of preferred fiber directions and may thus 
have interesting underpinnings in biomechanical applications. Therefore it is 
possible to learn a lot from a simple semi-inverse problem about the complex 
structure of nonlinear elasticity. 

We conclude with an open problem. For the incompressible Mooney- 
Rivlin material, Rajagopal and Wineman (1985) found solutions to the bal- 
ance equations within the class described by f fT7|) . with ffl8|) in force when 
Q = ipZ + ip , i.e. when the angles of rotation of the two bounding planes of 
the slab are different. We ask whether it is possible to find a special class of 
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compressible materials for which this is possible? 
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